NONLINEAR DYNAMO IN A SHORT TAYLOR-COUETTE SETUP 



C. NOREl•2'^ J.-L. GUERMOND*'*, R. LAGUERRE^'^ J. LEORAT^ F. LUDDENS^ 

Abstract. It is numerically demonstrated by means of a magnetohydrodynamics code that a short 
Taylor-Gouette setup with a body force can sustain dynamo action. The magnetic threshold is 
comparable to what is usually obtained in spherical geometries. The linear dynamo is characterized 
by a rotating equatorial dipole. The nonlinear regime is characterized by fluctuating kinetic and 
magnetic energies and a tilted dipole whose axial component exhibits aperiodic reversals during 
the time evolution. These numerical evidences of dynamo action in a short Taylor-Couette setup 
may be useful for developing an experimental device. 



1. Introduction 

Still a century after Larmor suggested that dynamo action can be a source of magnetic field in 
astrophysics, the exact mechanism by which a fluid dynamo can be put in action in astrophysical 
bodies remains an open challenge. In addition to the numerous analytical and numerical studies 
that have been done since Larmor's work, it is only recently that fluid dynamos have been pro- 
duced experimentally [U [131 [12]. These experimental dynamos have been helpful in particular to 
explore the nonlinear saturation regime. For instance, the dynamo produced in the Cadarache 
experiment [12j has an axial dipolar component and exhibits polarity reversals that are not unlike 
those observed in astronomical dynamos. The design of this experiment, however, has peculiar 
features that distinguishes it from natural dynamos. The most notable one is that the flow motion 
is induced by counter-rotating impellers. This driving mechanism induces an unrealistic differ- 
ential rotation in the equatorial plane and produces a large turbulent dissipation. Even with a 
mechanical power injection close to 300 kW, the magnetic Reynolds number of the flow of liquid 
sodium hardly reaches i?m = 45. Another peculiarity of this experiment is that dynamo action 
has not yet been obtained by using blades made of steel. The dynamo threshold has been reached 
at i?m = 32 by using blades made of soft iron instead. The objective of the present work is to 
investigate an alternative driving mechanism that shares the fundamental symmetry properties 
of natural dynamos, i.e., axisymmetry and equatorial symmetry (so-called SO(2)-Z2 symmetry). 
The Taylor-Couette geometry is a natural candidate for this purpose, since this configuration is 
already known to produce dynamos both in axially periodic geometries [15] and in finite vessels of 
large aspect ratio [7j. We examine in the present paper the dynamo capabilities of Taylor-Couette 
flows in vessels of small aspect ratio, and we compare the results obtained in this setting with 
those from more popular spherical dynamos [3]. 

The paper is organized as follows. The formulation of the problem and the physical setting 
of the Taylor-Couette configuration under consideration are described in Q The formulation of 
the problem and the physical setting of the Taylor-Couette configuration under consideration are 
described in ^ Three types of flows are considered in the paper and are discussed in ^ These 
flows are the standard Taylor-Couette flow driven by viscous stresses, a manufactured Taylor- 
Couette flow, and an optimized flow driven by a body force that models rotating blades attached 
to the lids. Two kinematic dynamo conflgurations are investigated in Q It is found that the 
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poloidal to toroidal ratio of the velocity field generated by viscous driving only (standard Taylor- 
Couette) is not large enough to generate a dynamo at Rm < 200. Dynamo action is obtained by 
using the strengthen Taylor- Couette flow and the forced Taylor-Couette flow. In both cases the 
poloidal to toroidal ratio of the velocity field is close to one. A nonlinear dynamo obtained with 
the forced Taylor-Couette flow is described in ^ In the early linear phase of the dynamo, the 
magnetic field at large distance is dominated by an equatorial rotating dipole. In the established 
nonlinear regime, an axial axisymmetric component of the magnetic dipole is excited and exhibits 
aperiodic reversals. Concluding remarks are reported in f|6j 

2. Formulation of the problem 

2.1. The physical setting. We consider an incompressible conducting fluid of constant density p 
and constant kinematic viscosity ly. This fluid is contained between two coaxial cylinders of height 
Lz- The radius of the inner cylinder is Ri and that of the outer one is Rq. The inner cylinder 
is composed of a solid conducting material. The inner cylindrical wall and the top and bottom 
lids corotate at angular velocity Oj. The outer cylindrical wall is motionless. The inner solid core 
may rotate or not, i.e., the inner core and the inner cylindrical wall may have different angular 
velocities. The conductivity of the fluid and inner solid is assumed to be constant and is denoted 
ctq. The magnetic permeability fiQ is assumed to be constant in the entire space. 

Let U he a, reference velocity scale yet to be defined. We then consider the following reference 
scales for length, C = Ro — Ri, magnetic field, H = U^/pJJjlq, and pressure, V = pU^ . The non- 
dimensional parameters of the system are the kinetic Reynolds number, i?e, the magnetic Reynolds 
number, i?m, the radius ratio rj, and the aspect ratio F: 

(2.1) Re = —, Rm = fJ'OcroUC, 11=-^^ r=^. 

To limit the number of geometrical parameters, we restrict ourselves in this paper to = 0.5 
and F = 2. Abusing the notation, this immediately implies that i?j = 1 and i?o = 2 in non- 
dimensional units. We did not explore other aspect ratios (see for example |141 [H [TTl [9] for short 
aspect ratios and different angular velocities). The conducting domain is partitioned into its 
fluid part enclosed between the two walls, ficj, and its solid part enclosed in the inner cylinder, 
Ocs. Using non-dimensional cylindrical coordinates {r,9,z), we have flcf = [Ij 2] x [0, 27r) x [— 1, 1] 
and 0,cs = [0, 1] x [0, 2tt) x [— 1, 1] . The conducting material is embedded in a non-conducting region 
denoted Q^, which we refer to as the vacuum region. 

The non-dimensional set of equations that we consider is written as follows in the conducting 
material: 

(2.2) dtu + {u-V)u + Vp = -^Au + (VxH")xff + f/ 

Rc 

(2.3) V-u = 

(2.4) dtW -VxiuxW) = — AH^ 

(2.5) VH" = 0, 

where u, p, and H*^ are the velocity field, pressure, and magnetic field, respectively. The magnetic 
field in il.^ is assumed to derive from a harmonic scalar potential: H'' = Vcj), A(j) = 0. The 
transmission conditions across the interface separating the conducting and nonconducting material 
are such that the tangent components of the magnetic and electric fields are continuous (see |6j). 

We consider three different settings: (i) The incompressible Navier-Stokes setting (H'^ = 0); 
(ii) The Maxwell or kinematic dynamo setting; (iii) The nonlinear magnetohydro dynamics setting 
(MHD). In the Navier-Stokes setting, H"^ is set to zero in the Lorentz force and the induction 
equation is not solved. The source term f/ is an ad hoc body force that models blades fixed at the 



endwalls, see ^3.3 When f/ = 0, the viscous stress induced by the rotating walls is the only source 
of momentum, see §3.1[ In the Maxwell setting, only the induction equation is solved assuming 
that some ad hoc velocity field u is given. In the MHD setting, the full set of equations is solved. 
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Since the definition of tlie reference velocity in similar dynamo configurations may be different 
(velocity at a given point, maximal speed in the flow, etc.), we introduce the root mean square 
(rms) velocity to facilitate comparisons: 

(2.6) W' = —^f ||u(x,t)fdx, 

vol(S2c/) J^^^. 

where the dimensionless fluid volume is vol(ilc/) = Svr in the present case. 

2.2. Numerical details. The code (SFEMaNS) that we have developed solves the coupled Navier- 
Stokes and Maxwell equations in the MHD limit in heterogeneous axisymmetric domains composed 
of conducting and nonconducting regions by using a mixed Fourier/Lagrange finite element tech- 
nique. Continuous Lagrange Finite elements are used in the meridian plane and Fourier modes 
are used in the azimuthal direction. Parallelization is done with respect to the Fourier modes. 
Continuity conditions across interfaces are enforced using an interior penalty technique [HI E]- 
SFEMaNS can account for discontinuous electrical conductivity and magnetic permeability distri- 
butions [5I, f8] . The magnetic field in the nonconducting regions is assumed to derive from a scalar 
magnetic potential, i.e., the configurations that we model are such that there is some mechanism 
that ensures that the circulation of the magnetic field along any path in the insulating medium is 
zero (this happens for instance when the vacuum is simply connected). The velocity field in rj^/ 
and the magnetic field in J7c are approximated using continuous P2 polynomials, and the pressure 
field in f^c/ is approximated using continuous Pi polynomials. In the vacuum Q,^, the magnetic 
potential (j) is approximated using continuous P2 polynomials. Typical characteristics of the meshes 
in the meridian section of all the cases studied in this paper are summarized in Table [T| 



Run 


Ax 


At 


np{P) 


np{V) 


np{H) 


np{4>) 


M 


d • 

^3.3 
ii.l 


5.1 

, \ 

5 


3.2 
4.2 


1/100 
1/100 
1/100 
1/100 


0.025 
0.025 
0.005 
0.005 


5911 
5911 
5911 
5911 


23341 
23341 
23341 
23341 


29821 
29821 


14041 
14041 


8 
12 
4 
32 



Table 1. Characteristics of the runs: Ax is the quasi- uniform meshsize in Qc', 
At is the timestep; np(P) is the number of Pi nodes for the pressure field in Oc/i 
np{V) is the number of P2 nodes for the velocity field in fie/; np{H) is the number 
of P2 nodes for the magnetic field in f^c; np{(f)) is the number of P2 nodes for the 
magnetic potential in Q^. The numbers np{P), np{V), np{H) refer only to the 
meridian section. The total number of grid points for each unknown Y is obtained 
by multiplying np(Y) by 2 times the number of Fourier modes, M, minus one. 



The performance of SFEMaNS has been validated on various kinematic and nonlinear dynamo 
configurations. In particular, a study of two Taylor-Couette setups using SFEMaNS is reported 
in [7J. In the first case F = 4, 77 = 0.5, and z-periodicity is assumed; in the second case F = 
2tt, rj = 0.5 and the vessel is finite, i.e., no z-periodicity is assumed and the vessel is closed at 
both ends. In both cases the inner wall rotates, but the outer wall and the two endwalls (when 
present) are motionless. The self-consistent saturated dynamo found in [15j in the z-periodic case 
has been reproduced in and a new nonlinear dynamo has been found in the finite vessel at 
Re = 120, i?m = 240. The behaviors of the z-periodic and finite-vessel dynamos, as observed in 
[7], significantly differ. After some transient, the kinetic and magnetic energies of the z-periodic 
dynamo converge to a stationary value. The final nonlinear MHD state is a steady rotating wave 
resulting from the balance between the driving effect of the viscous shear and the braking effect 
of the Lorentz force. The nonlinear dynamo action found in the finite vessel shows a different 
behavior in which the spatial symmetry about the equatorial plane (or mid-plane) of the velocity 
and magnetic fields plays a key role. The dynamo is cyclic in time and the fields rotate rigidly 
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with modulated amphtude. In these two cases (periodic and finite extension), the wavelength 
of the magnetic eigenvector is about twice that of the flow; as a result, the velocity field in the 
median plane of a single magnetic structure is directed inwards. This feature is shared by the 
spherical kinematic dynamos studied in [3]. It is reported in [3] that the lowest critical magnetic 
Reynolds number is obtained when the velocity field forms two poloidal cells that flow inwards in 
the equatorial plane. Note in passing that the two numerical experiments reported in [7] clearly 
confirm that assuming periodicity or enforcing finite boundary conditions give rise to dynamos with 
fundamentally different behaviors, i.e., assuming periodicity or ad hoc boundary conditions for the 
sake of numerical convenience may have nontrivial consequences. The series of observations above 
have led us to investigate more thoroughly the Taylor-Couette configuration with aspect ratio 

r = 2. 

3. Hydrodynamic forcing 

Since a number of dynamo studies have shown that the ratio of poloidal to toroidal speed should 
be close to unity to obtain the lowest critical magnetic Reynolds number, it is important to control 
this ratio. We describe in this section the mechanisms that we use to optimize the velocity field 
for dynamo action. 

3.1. Taylor-Couette flow (viscous driving only). When the aspect ratio is about 2 and the 
kinematic Reynolds number is moderate, two counter-rotating poloidal cells form with a toroidal 
angular velocity oriented in the same direction as that of the inner cylinder. In order to enforce 
the equatorial jet to flow inwards, we let the lids of the vessel rotate with the angular velocity of 
the inner cylinder and we keep the outer cylinder motionless. Note that it is important to have 
the lids and the inner cylindrical wall of the vessel to corotate; this makes the equatorial jet flow 
inwards and makes the overall velocity filed similar to the spherical flows that are known to yield 
dynamo action [3]. 

We define the velocity reference scale to be 

(3.1) U = VliRi, 

when the only source of momentum is the viscous stress at the boundary. 




Figure 1. Taylor-Couette flow Vo, F = 2, iic = 120. The angular velocity of 
the lids and the inner cylinder is fi, = 1; the outer cylinder is motionless. A 
radial jet flows inward at the equator. Represented are the poloidal flow (vectors 
and streamlines) —0.2 < Vr < 0.3, and —0.25 < Vz < 0.25, and the toroidal (or 
azimuthal) flow 0.25 <V0<2 (every 0.25). 

At i?e = 120 in the Navier-Stokes regime, the flow is steady, and forms the expected two toroidal 
cells invariant under the SO(2)-Z2 symmetry, i.e., axisymmetric and symmetric with respect to the 
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equatorial plane, see Figure [T| This flow, henceforth generically referred to as Vq, is characterized 
by its rms velocity, Vq, defined as follows: 

where and Vq^ are the rms poloidal and toroidal velocities of the reference hydrodynamic flow, 
respectively, and vol(r2c/) = Gvr is the volume of the vessel. Our computations give Vq = 0.272; 
this value is significantly lower than the maximum speed at the rim of the endwalls which is equal 
to 2. The poloidal to toroidal ratio is Aq = V^op/^l ~ 0.235. 

We have verified that the flow Vq is stable with respect to non-axisymmetric perturbations 
supported on the Fourier modes m = 1, • • • , 7 at i?e = 120. The velocity field Vq together with a 
sketch of the setup is shown in Figure [l} This reference flow is used in §4.1| to perform kinematic 
dynamo simulations . 

3.2. A modified Taylor-Couette flow. In order to perform kinematic dynamo simulations with 
a velocity held that has a poloidal to toroidal ratio that can be controlled easily, we construct an 
ad hoc field based on Vq. We use the poloidal and toroidal components of the vector field Vq to 
define a kinematic field, V^, with a pre-assigned poloidal to toroidal ratio as follows: 

(3.2) V, = -^^op + ^rVot. 

a[e) a{e) 

The normalization is done so that the rms of is the same as that of Vq. This gives 

(3.3) a\e)=^-^^, A{e) = eAo- 

i -I- ii-o 

Since the toroidal component of the velocity at the inner cylinder is equal to l/a(e), the angular 
velocity of the inner wall is fij = l/a{e), and this also means that the reference velocity scale is 

(3.4) U = a{e)QiRi. 

Although the vector field is not a solution of the Navier-Stokes equations, it is nevertheless 
solenoidal. This fiow is henceforth generically referred to as Ve- Computations have been done 
(see ^.1) for the values of e reported in Table [2] The quantity denoted Fmax in Table [2] is the 



maximum of the velocity modulus; Vmax depends on e. 
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0.71 


0.94 


1.18 


1.41 
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1.67 


1.49 


1.32 


1.20 
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1.23 


1.25 
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Table 2. Modified Taylor-Couette fiow: normalization factor a(e), poloidal to 
toroidal ratio A(e) and maximum of the velocity modulus 14nax- 



3.3. Forced Taylor-Couette flow (viscous driving plus body force). A number of dynamo 
studies have shown that the ratio of poloidal to toroidal speed should be close to unity to obtain 
a low critical magnetic Reynolds number. Viscous driving by the rotating walls yields a value for 
this ratio that is not close to unity (Aq = 0.235 at R^. = 120, see section above). At low Reynolds 
numbers, the fiow is steady and axisymmetric. It is relatively easy to vary the relative amplitude 
of the toroidal component in experimental setups by using blades fixed to the corotating endwalls 
to act as centrifugal pumps. This configuration, however, is difficult to implement in a computer 
code. In order to better control the poloidal to toroidal ratio in our simulations, we have chosen 
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to model the toroidal driving by a body force. The action of blades on the top and bottom lids is 
modeled by an ad hoc axisymmetric divergence- less force given in dimensional form as follows: 

A 

P- 



(3.5) 



< 1 and 1.2 < r < 1.^ 



if 0.8 < |z, 
otherwise, 

where the non-dimensional parameter A has been tuned to optimize the poloidal to toroidal ratio. 



Note that (3.5) defines the reference velocity U. The resulting velocity field is denoted V/ and the 
flow is generically called V/. 
We have found that using A 



2.5 at Rc = 120 gives Aj 



1.04 and the rms velocity 



is Vj = 0.219. We have observed that the azimuthal velocity in the vicinity of the inner radius 
is close to 0.55; hence, to reduce the viscous boundary layer at the inner wall and endwalls, we 
have set the dimensionless angular velocity to Vli = 0.55. The steady axisymmetric flow Vi is 
shown in Figure [2| (see Figure [T] for a comparison with the pure Taylor-Couette flow) . We have 





(a) Poloidal flow, 9 = n 



(b) Toroidal flow, 9^0 



Figure 2. Forced Taylor-Couette flow Vi, T = 2, Rc = 120, and A = 2.5. The hds 
and the inner cylinder rotate with angular speed = 0.55; the outer cylinder is 
motionless. A radial jet flows inwards at the equator. Represented are the poloidal 
flow (vectors and streamlines), —0.7 <Vr< 1.4, —1.1 <Vz< 1.1, and the toroidal 
(or azimuthal) flow, < < 1.6 (10 levels) 

verified, by performing nonlinear Navier-Stokes simulations, that the flow V/, at Re = 120, is stable 
with respect to three-dimensional perturbations supported on Fourier modes up to m = 11. The 
first hydrodynamic non- axisymmetric instability occurs on the Fourier mode m = 3 at Re = 168. 



The steady and axisymmetric forced Taylor-Couette flow Vi is used in ^4.2 to perform kinematic 
dynamo simulations. 

3.4. Summary. To compare the flows Vq, Vi, and Ve, we show in Table[3]the following character- 
istics of these three flows: rms velocity, V*; maximum of the velocity modulus in the fluid domain, 
Vmax! poloidal to toroidal ratio, A. 



Run 


V* 




A 


Viscous flow Vo, 
Modified flow Ve, V- 
Forced flow Vi, ^3.3 


^3.1 

5 2,^ 
,U.2 


o 

]5| 


0.272 
0.272 
0.219 


2.00 
(1.20,2.00) 
1.09 


0.235 
0.235x6 
1.04 



Table 3. Characteristics of the flows Vo, Ve and Vj: V* is the r.m.s. speed; Fmax 
is the maximum of the velocity modulus in the fluid domain; A is the poloidal to 
toroidal ratio. 
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4. Kinematic Dynamos 

We evaluate in this section the properties of the kinematic dynamos generated by the flows Ve 
(viscous driving) and V/ (viscous driving plus body force). 

4.1. Parametric study of the poloidal to toroidal ratio using Ve- We investigate the dynamo 



properties of the manufactured flow Ve in the kinematic regime, see ^3.2 The reference velocity 



scale is defined in (3.4). To ensure that the velocity is continuous across the solid/fluid interface, 
the angular velocity of the inner core is set to be l/a(e). The conductivities of the solid inner core 
and the fiuid are identical. 

We perform two studies at i?m = 100 and Rm = 200 to determine the optimal weight e that 
gives the largest growthrate of the dynamo action in the kinematic regime. The computations are 
done with SFEMaNS in Maxwell mode. The magnetic field is initialized to some small random 
values and the growth rate (i.e., the real part of the leading eigenvalue) is computed by running 
short time simulations for various ratios e G [3, 16] shown in Table [2j As the vector field is 
axisymmetric, the term V x (V eXH^) cannot transfer energy between the Fourier modes of H*^, 
i.e., the Fourier modes are uncoupled. The first bifurcation is of Hopf type and the most unstable 
eigenvector is the Fourier mode m = 1. The growthrate of the magnetic energy is reported in 
Figure [3j There is no dynamo action at Rm = 100. Dynamo action occurs at Rm = 200 in the 
range 4.2 < e < 15.4, which corresponds to 1.0 < A(e) < 3.8. Note that the purely viscous driving, 
which corresponds to e = 1 and Aq = 0.235, cannot sustain a dynamo at Rm = 100 and Rm = 200. 




-0.05 I ' ' ' ' ' ' 1 

2 4 6 8 10 12 14 16 



e 

Figure 3. Kinematic dynamo. Growthrate of Fourier mode m = 1 for the modified 
Taylor-Couette fiow, V^, as a function of e for R^ = 100 and i?m = 200; T = 2 and 
Re = 120. 

Figure |4] shows the magnetic eigenvector for the Fourier mode m = 1 at R^ = 120, i?m = 200 
and eopt = 8. This eigenvector is a rigid wave that rotates in the same direction as the inner 
cylinder and top/bottom lids, and its period of rotation is T ~ 870, i.e., more than 66 rotation 
periods. Upon introducing the equatorial symmetry operator 5z2H = {Hr, Hq, —Hz){r,6,—z). 
the magnetic field has the following symmetry property: 

(4.1) = 5z2H^ 

i.e., the magnetic field has the same symmetry as the velocity field. 

We now evaluate the critical magnetic Reynolds number and its minimal value with respect to 
e. We assume that the growthrate depends smoothly on R^. Upon inspecting Figure [3] we see that 
the growth rate is maximum for Copt ~ 5 at i?m = 100 and for eopt ~ 8 ate Rm = 200. Then, by 
drawing the line connecting the two maximum points in Figure [3| we observe that this line crosses 
the horizontal line of zero growth rate in the interval e G [6.5, 7]. We have chosen to explore the 
value e = 6.5, which gives the poloidal to toroidal ratio A = 1.53. The growth rate for various 
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(d) Hr aie=l (e) He &t 9 = ^ (f) a.t 6 ^ I 



Figure 4. Kinematic dynamo with flow V(e = 8) at Re = 120, i?m = 200. Mag- 
netic eigenvector for Fourier mode m = 1. Represented in (a) to (f) are the radial, 
azimuthal, and vertical components, normahzed by the square root of the magnetic 
energy, in two complementary planes, with 0<r<2, — 1<2;<1 (the z— axis 
is on the left): for 6* = 0, -0.85 < Hr < 0.85 (every 0.17), -0.1 < He < 0.68 
(every 0.17) and -0.85 < H - z < 0.85 (every 0.17); for 6 = ^, -1 < Hr < 3.75 
(every 0.25), -1 < He < 0.15 (every 0.25) and -1.5 < H^ < 1.5 (every 0.25). 
Represented in (g) is the isosurface ||H^|p (7% of maximum value) colored by the 
azimuthal component. Note the m = 1 structure. 



0.02 




-0.01 I ' ' ' ' ' ' ' 1 

130 140 150 160 170 180 190 200 



Figure 5. Kinematic dynamo of flow V(e = 6.5), F = 2, i?o = 120. Growthrate of 
the Fourier mode m = 1 as a function of i?ni- 



magnetic Reynolds numbers in the range [125, 200] has been computed. The results are shown in 
Figure [5j We have found that the critical magnetic Reynolds number for the Fourier mode m = 1 
is Rtn = 155 with e = 6.5 and A = 1.53. 



4.2. Kinematic dynamo in the forced Taylor- Couette setup. We use the steady axisym- 
metric forced flow, Vi, at = 120 to perform kinematic dynamo computations. The reference 
velocity scale is defined in (3.5) with A = 2.5. To determine whether the rotation of the inner 
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solid core has any impact on the dynamo threshold, we have compared growth rates when the 
solid inner core is motionless and when the inner wall and solid inner core corotate with angular 
speed fij. Whether the inner core rotates or not, the magnetic eigenvector m = 1 is always the 
most unstable. We show in Figure [6] the computed growth rates in the two cases in the range 
i?m G [175,200]. The motion of the inner core does not seem to have a significant influence; in 
both cases the dynamo threshold is Rm = 180 for the Fourier mode m = 1. 
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Figure 6. Kinematic dynamo with V/ flow, F = 2 and Rc = 120 . Growthrate of 
the Fourier mode m = 1 as a function of i?ni- ROT: rotating inner core; NO-ROT: 
non-rotating inner core (but inner wall rotates). 



The structure of the magnetic eigenvector corresponding to the Fourier mode m = 1 with a 
rotating inner core is shown in Figure [7] (at Rg = 120, R^ = 200). It is a rigid wave that rotates 
in the same direction as the inner cylinder and top/bottom lids with period T ~ 120. 

5. Nonlinear dynamo action 

We report in this section on nonlinear dynamo computations done with the forced Taylor- 
Couette setup at i?e = 120 and R^ = 200. 

5.1. Description of the setting. We consider the forced Taylor-Couette setup described in 
^4.2 The reference velocity scale is defined in ( |3.5[ ) with A = 2.5. We perform nonlinear dynamo 
computations with the parameters Re = 120, i?m = 200. The inner core is kept motionless. We 
work with 32 azimuthal modes (m = 0, . . . ,31), and the meridional finite element mesh is the 
same as in the kinematic runs. The total number of degrees of freedom to be updated at each time 
step is 11,353,104. The initial velocity field is the axisymmetric flow Vi that we computed in the 
Navier-Stokes regime at -Re = 120. The initial magnetic seed is the growing Fourier mode m = 1 

200 > R^c 



180. 



obtained in the kinematic computations described in ^4.2 at R^ 

When dynamo action occurs, the magnetic energy grows exponentially until the Lorentz force 
is capable of modifying the base flow. This transient phase lasts about 5 rotation periods. When 
the Lorentz force is strong enough, a new regime settles where the magnetic energy saturates. 
Nonlinear saturation is a slow process that lasts at least 200 rotation periods (see Figure 11 in [7]). 
Although SFEMaNS is parallel with respect to the Fourier modes, the volume of computation 
required by this type of simulation is such that we have not been able to explore other kinematic 
and magnetic Reynolds numbers within the resources allocated to this project. The nonlinear run 
presented in this section used about 15600 cumulated CPU hours with 32 processors on an IBM 
Power 6 cluster. 

5.2. Time evolution of the energy. The time evolutions of the kinetic and magnetic energies 
are reported in Figure [8][^a-b), where the kinetic and magnetic energies are defined as follows: 
h In f II'^IP*^^! 5 In l|H^|pdx, respectively. From t = to t = 500 (first transition), the kinetic 
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(d) Hrak.e=l (e) He at 6 ^ ^ (f) iJ, at 61 = f 



Figure 7. Kinematic dynamo with flow V/ at Re = 120, i?m = 200. Magnetic 
eigenvector for Fourier mode m = 1. Represented in (a) to (f) are the radial, 
azimuthal, and vertical components, normalized by the square root of the magnetic 
energy, in two complementary planes: for 9 = 0, —0.9 < Hr < 0.2 (every 0.1), 
-1.4 < He < 0.35 (every 0.25) and -0.6 < < 0.6 (every 0.1); for = 7r/2, 
-2.2 <Hr< 0.9 (every 0.25), -0.25 < Hg < 1.75 (every 0.25) and -1.25 < < 
1.25 (every 0.1). Represented in (g) is the isosurface HH'^lp (14% of maximum 
value) colored by the azimuthal component. Note the m = 1 structure. 

energy decreases and the magnetic energy grows exponentially with a growthrate similar to that 
of the kinematic dynamo. Then both the magnetic and the kinetic energies seem to saturate in a 
first nonlinear regime, 500 ^ t < 1100. During the first transition, the fiuid flow loses the axial 
symmetry, m = 0, and the magnetic field loses the symmetry associated with the Fourier mode 
m = 1. The flow being forced by the Lorentz force (VxH^)xH'^, the velocity thereby acquires 
a contribution on the Fourier mode m = 2. The magnetic field being deformed by the action 
of the induction term uxH'^ acquires a contribution on the Fourier mode m = 3. The cascade 
of nonlinear couplings generate even velocity modes, m = 0,2,..., and odd magnetic modes, 
m = 1,3, .. . During this transitional phase that consists of populating the Fourier modes, the 
axisymmetry of the velocity field is broken but the equatorial (mid-plane) symmetry is preserved 
for both the velocity and the magnetic field. 

A second nonlinear transition starts at t = 1100 and lasts until t = 1175. In this time interval 
the magnetic energy increases and the kinetic energy decreases. This change of behavior is due to 
the breaking of the equatorial symmetry. This phenomenon is well illustrated by computing the 
energy of the symmetric part, ^{H'^ + 5z2H'^), and anti-symmetric part, ^(H.'^ — 5z2H'^), of the 
magnetic field. The time evolution of these two quantities is shown in Figure [8][^c). The equatorial 
symmetry breaking is driven by the small even azimuthal modes of the magnetic field as can be 
seen on Figure [9] (a)(b), especially the magnetic mode m = 2. 

In the time interval 1175 < t < 1600, the system enters a second nonlinear regime characterized 
by large fluctuations and a dynamics dominated by the large Fourier modes. Between t = 1600 
and t = 2200, we observe a third transition during which the small even modes of the magnetic 
field increase again until they reach the final saturated state. A third and final nonlinear regime 
settles beyond t = 2200. Figure loFc) shows that the large Fourier modes, exemplified by m = 30 
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Figure 8. Nonlinear dynamo in the forced Taylor-Couette setup, (a) Time evolu- 
tion of kinetic (-0.37) and magnetic energies in the conducting region < r < Rq 
and — r/2 < z < T/2. Panel (b) shows the symmetric and anti-symmetric compo- 
nents of the magnetic energy. 
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Figure 9. Time evolution of different modal energies; (a) magnetic energies m = 
0, 2 (bottom curves) and m = 1, 3 (top curves); (b) kinetic energies m = 1, 3 and 
magnetic energies m = 0, 2; (c) magnetic energies m = 1 (top curve), m = 30, 31 
(middle curves) and m = 0. 



and m = 31, basically fluctuate within some asymptotic range for t > 1250, whereas the small 
even modes grow until they become energetically significant. 



We show in Figure 10(a) the time series of the azimuthal component of the magnetic field at 
the point (r = 1.2, 6 = 0,z = —0.5). The envelop of the signal first grows then reaches a maximum 
range. The period of the signal is T ~ 112.5; this corresponds to a wave that rotates in the same 
direction as the inner cylinder and top/bottom lids. More frequencies appear beyond t = 1100; the 
signal is the superposition of an oscillation of period T f« 150 and a modulation of period Tmod ~ 17, 
which happens to be of the same order as the wall rotation period Tuds = 27r/0.55 ~ 11.4. The 
breaking of the equatorial symmetry is measured by monitoring the anti-symmetric part of the 
magnetic field. We show in Figure |10(b)| the time evolution of twice the anti-symmetric part of 
Hg at (r = 1.2,6* = 0,2; = -0.5). 

5.3. Spatial structure of the dynamo. To have a better understanding of the structure of 
the velocity and magnetic fields during the three nonlinear regimes identified above, we show 



in Figure 11 the isosurfaces of the magnetic and kinetic energies at t = 1000, t = 1400, and 



t = 2500. Since the axisymmetric velocity mode m = is dominant (see Figure [Tl[d-f)), we show 
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Figure 10. Time series of (a) Hg at (r = 1.2,6' = 0,2; : 
anti-symmetric part of Hq at (r = 1.2, = 0, z = —0.5). 



-0.5) and (b) twice the 



in Figure [TT](g-i) the kinetic energy without its axisymmetric contribution to better distinguish 
the fine structures. The magnetic field at t = 1000 is dominated by the odd azimuthal modes (see 
Figure [TT| a)) and resembles the eigenvector shown in Figure 7(g) The velocity field is composed of 
even modes as can be seen on Figure [TT|g) , see e.g. the two dark structures that are diametrically 
opposed. Similar spatial distributions are observed at time t = 1400 with the addition of smaller 
scales. At time t = 2500 the small scale modes are even more apparent, and the magnetic and 
velocity structures in the fluid domain are more deformed. 

To characterize the long distance influence of the magnetic field, we have recorded the time 
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shows the time 



evolution of the magnetic dipole defined by D = rx(VxH'^)dx. Figure 
series of the three Cartesian components of the magnetic dipole in the time interval < t < 2500. 
During the first two transitions and nonlinear regimes, i.e., < t < 1600, the dipolar moment 
is purely equatorial and rotates at the same frequency as the magnetic field. The axial moment 
starts to grow at the beginning of the third transition {t > 1600) and changes sign several times 
afterward (note that the time series in (a) is under sampled in this range). Figure 12 b) presents 
a zoom of the time evolution of Dz showing two reversals. The magnitude of the quadrupolar 
moments (data not shown) stay below 10~^ until t = 1175, then increase and saturate to values 
four times smaller than the magnitude of the dipolar moment. Magnetic field lines in the vacuum 
at t 



2500 show a pattern characteristics of an equatorial dipole (see Figure 13). 



6. Concluding remarks 

We have numerically demonstrated that pure viscous driving by smooth rotating walls in a 
short Taylor-Couette setup does not lead to dynamo action for < 200 since the poloidal to 
toroidal ratio of the velocity field is too small. An adjustment of the poloidal to toroidal ratio is 
needed to achieve dynamo action in the kinematic regime. We have implemented an ad hoc body 
force to produce a poloidal to toroidal ratio that is of the same order as what is needed in the 
kinematic simulations to trigger the dynamo action. This force may be thought of as a model for 
the action of blades fixed to the static or to the rotating walls/lids at convenient angles. This 
force has also the same symmetry properties as the geodynamo, i.e., the SO(2)xZ2 symmetry 
(axisymmetry and equatorial symmetry). The critical magnetic Reynolds number of this setup 
based on the inner cylinder speed fiQa^iRj = 0.55x180 = 99 is in the range of what has been 
obtained in the kinematic studies of [3J in a spherical container with the same symmetries. This 
magnetic Reynolds number is also comparable to what has been reported in [7j for Taylor-Couette 
simulations in vessels of larger aspect ratios and with pure viscous driving. 

A nonlinear simulation has been performed at Rg = 120, -Rm = 200 over 225 rotation periods. 
In the early linear phase of the dynamo, the external field is dominated by an equatorial rotating 
dipole. In the established nonlinear regime, an axial axisymmetric component is excited and 
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Figure 11. Nonlinear dynamo at t = 1000, t = 1400 and t = 2500, for Re = 120, 
Rra = 200: (a-c) isosurface of HH'^p (25% of maximum value); (d-f) isosurface of 
||V|p (25% of maximum value); (g-i) isosurface of ||V|p without the axisymmetric 
mode (10% of maximum value); color scale proportional to azimuthal component. 



exhibits reversals. The relation between the main flow parameters of the time-dependent angle 
formed by the dipole and the rotation axis calls for further investigations, since it is a basic feature 
of observed planetary dynamos. 

Exploring the feasibility of an experimental fluid dynamo based on the present design will re- 
quire expertise from many different experimental and numerical fields [10]. To achieve a magnetic 
Reynolds number equal to 100 in a flow of liquid sodium requires that the kinematic Reynolds 
number be of order 10^. It is well known that such a value corresponds to a highly turbulent flow 
that can be studied only in experimental facilities, since it is far beyond the capacity of direct 
numerical simulations. The objective of such experiments should be to recover optimized poloidal 
and toroidal components after time averaging, which presumably would guide the design of the 
blades fixed to the endwalls. These experiments would also inform about the power requirements. 
Using a standard rotation frequency of 50 Hz, a magnetic Reynolds number of 100 can be ob- 
tained in liquid sodium at 150°C with an inner radius of approximatively 18 cm and an outer 
radius and height of 36 cm. This seems feasible since these dimensions are not far from those of 
the Cadarache experiment [12]. We conjecture however that the power required by this experiment 
at a given rotation frequency should be smaller, since the turbulence rate induced by co-rotating 
lids/impellers should be smaller than that of counter-rotating lids/impellers. A dynamo facility 
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(a) < t < 2500 (b) Zoom in 2410 < t < 2570 

Figure 12. Time evolution of the dipolar magnetic moment. Note that the vertical 
component has been multiplied by 5 in the left panel. 




(a) From above (b) From a perspective point of view 



Figure 13. Isosurface of ||H^|p (25% of maximum value) and magnetic field lines 
at t = 2500. 

presenting similarities with the present proposal is currently investigated by Colgate and collab- 
orators t2j. Their MHD device uses also a Taylor-Couette forcing in a short cylindrical container 
with size and targeted magnetic Reynolds number similar to those studied in the present paper. 
There are however differences: the flow in their experiment forms an outwards jet in the equatorial 
plane and is driven by viscous stresses only. More detailed comparisons of the respective merits of 
both designs should certainly be instructive. 
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